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Plastic  Soil  Mechanics  Theories. 


1.  Introduction.  Although  civil  engineers  have  been  inter- 
ested in  soil  mechanics  for  ages,  only  recently  has  It  begun  to 
be  feasible  to  apply  to  it  the  techniques  of  modern  plasticity 
theories.  There  is  little  experimental  evidence  as  to  which  theory 
to  choose.  On  the  other  hand  there  is  still  relatively  little 
theoretical  guidance  as  to  fruitful  lines  of  experimental  attack. 
This  report  is  an  attempt  to  bring  the  theoretical  and  experimental 
forces  somewhat  nearer  common  ground. 

To  be  mathematically  tractable  a theory  of  soil  mechanics  must 
drastically  simplify  true  soil  structure.  Only  a few  important 
soil  mechanisms  may  be  considered.  We  will  adopt  rather  naive 
models  of  these  mechanisms,  for  simplicity  is  at  a premium  and 
high  accuracy  is  not. 

The  basic  soil  structure  is  its  skeleton  of  solid  sand  or  silt 
grains,  usually  of  quartz.  They  are  elastic  in  nature.  Although 
they  are  crystalline,  their  crystallographic  axes  are  usually  ran- 
domly oriented  so  that  the  skeleton  as  a whole  may  be  regarded  as 
isotropic  .«■ 

The  Interstices  between  the  particles  are  filled  with  water, 
air,  or  both.  The  water  contains  flocculated  colloidal  particles 
much  smaller  than  the  sand  or  silt  grains,  but  appearing  more  pro- 
minently in  clays  than  in  sands.  They  are  subject  to  electro-chemi- 
cal forces  tending  to  bond  together  the  grains  of  the  skeleton. 

# In  some  alluvial  deposits  stream  action  has  caused  the  deposit 
of  non-randomly  oriented  highly  flattened  particles.  The  present 
theory  does  not  apply  to  such  cases . 
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A similar  bonding  is  caused  by  surface  tension  in  the  water  when 
air  is  present.  We  will  account  for  this  bonding  by  introducing 
a cohesive  pressure  pc , constant  or  at  most  weakly  dependent  on 
the  plastic  work  W,  The  energy  loss  resulting  from  the  breaking 
of  these  bonds  and  the  rubbing  together  of  the  skeletal  particles 
as  the  soil  deforms  will  be  accounted  for  by  a mechanism  of 
Coulomb  friction. 

The  soil  skeleton  will  transmit  a stress  . In  addition 

the  water  will  transmit  a pore  pressure  u,  dependent  primarily 
on  the  strain  6.  In  many  cases  the  cohesive  pressure  and  the 
pore  pressure  can  be  lumped  together  into  the  intrinsic  pressure 

P = u - pc. 

A soil  may  be  considered  incompressible  only  if  it  is  saturated 
with  water  and  contains  no  air  and  if  the  water  is  not  free  to 
flow  out. 

In  the  first  13  sections  several  different  plastic  soil 
mechanics  theories  are  developed.  In  section  14,  the  general 
question  of  testing  these  theories  is  discussed,  and  specialized 
in  section  15  to  the  ordinary  triaxial  type  of  testing  situation. 

In  section  16  the  theory  of  section  15  is  extended  for  the  perfect 
plasticity  case,  and  a numerical  example  is  given  in  section  17, 

Rather  accurate  experimental  work  is  required  in  testing  the 
theories,  as  one  might  expect.  If  the  experiments  suggested  are 
feasible  they  should  definitely  throw  light  upon  the  internal 
mechanism  of  soil  plasticity. 
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2.  The  total  stress  tensor  CT^j  and  the  intrinsic  pressure 
p are  defined  in  terms  of  the  skeletal  stress  tensor  tiie 

cohesive  pressure  p , and  the  pore  pressure  u by 

<ru  = Tlj  - P 5lj>  p = “ - pc  (1) 

where  is  the  Kronecker  delta:  5^  = 1,  5 ^ j = 0 for 

i t J. 

The  usual  mechanical  equations  of  equilibrium  of  plasticity 
and  elasticity  are 

* * Pl  = 0 <2> 

a x* 


where  the  summation  convention  is  used,  and  where  is  the 
rectangular  x^-component  of  external  force  per  unit  volume. 
The  total  strain  tensor  is  defined  by 

_ X ( 3ui  + ) 

eij  £ \ 5xj  dx^  J * 


(3) 


where  u^  is  the  component  of  displacement  in  x^-direction. 


This  may  be  resolved  into  an  elastic  strain  tensor  G®j  and 
plastic  strain  tensor  G^j  such  that 


hi  - eh  + e?j 


(4) 


By  [1;  §§  67,  69]*  the  elastic  part  is  given  by 


d6ij  = X d *ij  “ fi  5 ijd  ^~kk» 


(5) 


* Numbers  in  brackets  refer  to  the  references  at  the  end  of  the 
report . 
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where 


ec  = 1/(2 fJL)  = (1  + V )/E,  /3  = \/{2fJL^K  + 2/Ll)]  = 1//E,  (6) 

where  A ,yU.  are  Lame  constants,  E is  Young’s  modulus,  and  U 
is  Poisson’s  ratio. 

3.  Yield  conditions.  Yield  conditions  for  soils  have  usually 
been  based  on  the  notion  of  Coulomb  friction,  as  mentioned  in 
§ 2,  beginning  with  the  work  of  Coulomb  himself  [2],  and  have 
been  very  extensively  tested.  At  any  given  point  in  a soil  the 
skeletal  stresses  acting  on  an  arbitrary  plane  through  that 
point  may  be  obtained  by  transformation  from  the  stress  components 
referred  to  fixed  axes.  The  ratio  of  shearing  stress  to  normal 
stress  is  maximized  by  varying  the  direction  of  this  plane. 
According  to  the  Coulomb  friction  idea,  this  maximum  must  remain 
below  a certain  constant,  tan^  , ^ being  the  angle  of  shear- 

ing resistance.  It  is  a straightforward  calculus  problem  to  prove 


that 


1-sln  4 
1+s  in  j , 


(7) 


where  T^,  “TJi*  , are  the  three  principal  skeletal  stresses, 

so  arranged  that 

Tm<  % < T„  < 0 (8) 

the  convention  being  that  compressive  stresses  are  negative. 

(7)  may  be  written  in  the  form 

s < c + (T^tan , (9) 

where 


• = ~ rm)  • r'  ='pc  - T»  * Tm)  * 0 = pc  tan  (10> 
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(9)  is  the  standard  form  of  the  yield  condition  as  used  in  engin- 
eering soil  mechanics.  [3;  p.  22,  formula  [6]].  c is  called  the 
cohesion  of  the  material. 

The  yield  surface  corresponding  to  (7)  is  a six-sided  pyra- 
mid . 


(7)  is  simple  and  easy  to  apply  when  the  principal  directions 
are  known  from  symmetry  or  otherwise.  In  this  respect  it  is 
similar  to  the  Tresca  condition.  Like  the  Tresca  condition  it 
is  rather  complicated  when  expressed  in  invariant  form.  Then 
the  yield  condition  is 

xt(2x  + 3)  j3  - §(*+  i)  jxj2  + §r  * A)2 


(ID 


+ (3*  + 4)(j2  - ^j|)  = o, 


where 

h = rM  = T„  + + Tm, 

J2  = h rij  rij  “ - I sij  sij 

= F[(r«-  V2  + ( T«-  V)2  + ( V *n)2h  (12) 

J3  = l^ljl  + §J1J2  " ^7J1  = I Sljl 

= IJ1^  Tf*  “ 3Jl^Tnf  TJ1)» 

$ij  = Tij  ” 3J1  ®ij»  ^ = ftan2^  • 

Drucker  and  Prager  [4;  p.  153]  (see  also  [5;  p.  447])  have 
suggested  a simpler  yield  condition  equivalent  to 

jl/2  + 3-1/2^  jx  = 0 


(14) 
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for  constant.  This  has  a conical  yield  surface. 

The  homogeneity  of  conditions  (11)  and  (14)  makes  a two- 

dimensional  graphical  representation  possible.  Define 

x = _ 31/2j1/2/jij  y = (33/2/2)  J“3/2j3>  (15) 

Then  (11)  becomes 


X[(2 K + 3)y  + 3(X  + l)/x  - 

X/x3]2  + (3X  + 

4) (y2  - 1)  : 

which  may  be  solved  in  the  form 

-r  _ 3X(2X+3)  f -1 

Xx“3 

4(  x + 1 r L 

" 3 (K  + 1) 

+ + 4)  r 

i - **'2 ) //* 

K 1 

X2(  x + iyV2  L 

1 k+  ijy[ 

‘ 4 ( X + l)J 

(16) 


(14)  becomes 

X = (17) 

The  point  (x,y)  as  defined  by  (15)  can  cover  only  a 
finite  region  of  the  plane.  By  a straightforward  maximum- 
minimum  procedure,  which  will  be  omitted  here,  it  may  be  shown 
that 

x i 5 y ^ y i |x‘3  - I*"1. 

The  region  to  the  left  of  the  curves  (16)  or  (17)  is  elastic, 
and  the  region  to  the  right  is  forbidden.  During  plastic  deforma- 
tion the  point  (x,y)  lies  on  one  or  the  other  of  these  curves. 

In  Pig.  1,  equation  (16)  is  plotted  for  the  values  of  X indica- 
ted beside  each  curve  near  the  bottom.  These  curves  are  quite 
similar  to  those  for  (17)  for  small  K , corresponding  to  slight 
soil  compaction,  but  the  difference  becomes  progressively  more 
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marked  as  74  and  the  compaction  increase. 

4.  Strain  hardening.  The  angle  of  shearing  resistance 

is  fairly  constant  over  moderate  ranges  of  stresses  provided  the 
soil  is  somewhat  compacted  to  begin  with.  However  if  the  soil  i8 
initially  very  loose,  or  if  the  stress  becomes  high  enough  to 
drive  water  out  of  the  soil  permanently,  this  angle  may  change. 
This  phenomenon  is  referred  to  as  strain  hardening.  Following 
the  custom  in  metal  plasticity  we  will  regard  this  as  a function 
of  plastic  work  W,  where 

dW  = ^ijdeP  . (18) 

V/hen  strain  hardening  is  ignored  so  that  ^ is  considered  con- 
stant, the  resulting  theory  is  referred  to  as  perfect  plasticity. 

5.  The  pore  pressure.  If  possible  the  water  in  the  soil 
flows  In  such  a way  as  to  relieve  the  pore  pressure  at  the  expense 
of  the  skeleton.  However  the  inertia  and  viscosity  of  the  water 
oppose  this.  Similar  forces  act  on  the  air,  but  their  magnitude 
is  much  smaller  so  that  may  usually  be  neglected. 

By  the  law  of  Darcy, 

V OC  Vu, 

#5* 

^ being  the  velocity  of  flow  of  water  in  the  soil.  Now  V.* 
represents  the  outflow  of  water  per  unit  time  from  a unit  volume 
of  soil,  so  the  increase  per  unit  time  in  the  space  occupied  by 

g 

the  fluid  per  unit  volume  is  proportional  to  V u.  Two  factors 
contribute  to  this  increase  in  the  region  occupied  by  the  water. 
One  is  an  increase  dV  in  the  volume  of  the  trapped  air,  thus 


forcing  the  water  outward.  The  second  is  a decrease  equal  to 
-d®  = in  the  pore  space  within  the  skeleton.  This  decrease 

forces  the  water  outward  into  a larger  volume.  Therefore  we  have 

V2u  OC  Sv/3t  - Se/3t  (19) 

where 

e = (20) 

In  saturated  soil  V = 0.  In  unsaturated  soil  V > 0, 
but  may  be  taken  to  depend  on  the  pore  pressure  u.  Therefore 

h2  V2u  - 30/  d t = 0,  (21) 

where 

U = V(u)  - e,  v'(u)  < 0 (22) 

and  where  h is  essentially  constant,  but  may  tend  to  drop  off 
toward  zero  for  very  high  frequency  oscillations.  In  this  case 
it  might  be  considered  a function  of  % or  Vu, 

6.  The  plastic  strain  Increment . Drucker  [6]  has  shown 
that  there  exists  a plastic  potential  which  is,  in  fact,  the  yield 
function,  if  the  plastic  overwork  due  to  an  external  agency  acting 
on  a prestressed  material  is  always  positive.  While  this  is 
usually  the  case,  especially  in  metal  plasticity,  in  the  case  of 
the  yield  condition  (V)  it  may  not  be,  for  an  external  agency  may 
have  a "trigger”  action  in  which  it  releases  plastic  deformation 
while  gaining  energy  at  the  expense  of  the  initial  stressing  agency. 

We  will  assume  the  medium  to  be  isotropic,  which  seems  justi- 
fied if  the  soil  deformations  are  not  too  large.  We  will  also 
assume  the  plastic  strain  increments  dejj  to  be  proportional 
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to  the  plastic  work  increment  dW  defined  in  (18) , and  to  depend 
on  the  stress  components  an<^  possibly  also  on  u and  pc . 

Let  <T (Tg,  (T3  denote  the  principal  stress  components 
of  the  stress  tensor,  and  let 

Si  = tf'i  * 0*ec  (23) 


denote  the  principal  stress  deviations , Then,  letting  j = 
in  the  Appendix,  the  most  general  expression  for  the  plastic 
strain  increments  is 


deiJ  = (Eo  + E1  *ij  + E2  S1K  Skj>  dW> 

where 

Eo  = *2  *3  P1  + S3  si  P2  + *1  *2  p3> 

E1  = S1  P1  + % P2  - *3  P3 * <25> 

E2  = P1  * P2  + P3> 

where 

*1  = f(  <rlt  <rs,  cr3)/[(  sr  *8)(  *1-  *3)], 

Pg  = f(  (Tg,  (Tj,  cT^)/[(  Sg-  Sg-  s^)  ] , (26) 

FJ  = f(<r,,<r,,  <r,V[(  s,-  s,)(St-  s,)jt 

and  where  f(x,y,z)  is  arbitrary  except  that 

f(x,y,z)  = f(x,z,y),  (27) 

and  except  for  equation  (28)  below. 

If  we  multiply  (24)  by  (T ^ and  sum,  by  (18) 

(J1-3p)E0  + 2JgE1  + (3Jj  + fjjjg  . 2p Jg)Eg  = 1 
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Inserting  (25),  after  some  reduction, 

xf(x,y,z)  + yf (y,z,x)  + zf(z,x,y)  = 1.  (28) 

The  principal  axes  for  the  stress  are  also  principal  axes 
for  the  strain  increments.  Denote  principal  strain  increments  by 
dfiP,  deP,  d©P.  Then  by  (24)-(26), 

deP  = f ( <r1,  <r2,  <r3)dw,  de2p=  f(  <r2,  <r3,  <rx)dw. 


de3p=  f(  (r3^  <r1,  cr2)dw. 


(29) 


7.  The  Reus s equations . The  Reuss  equations  in  the  plastic 
theory  of  metals  are  derived  using  the  invariant  g as  a plastic 
potential.  For  these  equations 

f(  r1>  c r2,  <r3)  = ^(2^)  (30) 

8.  The  Tresca  theory.  When  the  Tresca  yield  condition  is 
used  in  conjunction  with  a plastic  potential, 

f*(  0“1,  <r2,  ars)  = [ 3gn(  0^-  0"2)  + sgn(  tf"1-  0“3)]/L(  - 0'2)sgn(^-^) 

(31) 

+ ( <r2  - <r3)  sgn(  <r2-  <r3)  + ( <r3-  d^sgn*  <r3-  o^)]. 


9,  The  Drucker-Prager  soil  theory.  In  this  theory  [4  ],  the 
yield  condition  (14)  is  used  in  conjunction  with  a plastic  potential. 
For  it, 

ftov  e-2,< r3)  = [3 B1/(84/2)I/[3’1/2'HU1-3p)  + #*,  (32) 
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10.  The  first  soil  theory.  In  this  paper  we  will  consider 
several  soil  theories.  In  the  first  of  these,  the  yield  condition 
(7)  is  used  in  conjunction  with  a plastic  potential.  This  leads 
to 

f(  0'1,  <r2 , <r3)  = [sgn(  - <r2)  + sgn(  er±  - orz)  * 2 sin^ 

- sin^  I sgn  (0^  - (T2)  - sg n ( CT’1  -^3)/ I/D, 

d = ((r1-<r2)sgn((r1-  <r2)  + ( r2-  <r3)sgn(  <r2-  <r3)  + 

, (33) 

+ ( 0"3  - CT^)  sgn(  ff”3  - o1^)  + 2 sin  p - 3p) 

- sin  <f>  [ <T1/sgn(  <T^  - er2)  - sgn  ( <T^  - 0^)/  + 0^/sgn(  (Tg-  cTj) 

- sgn(  CT 2~  0^)1  + <T3  / agn(  (T 3 - (7^)  - sgn(  ^ - 0g)/j. 

11.  The  second  soil  theory . If  one  assumes  that  the  soil 
shears  along  the  planes  of  maximum  ratio  of  shearing  stress  to 
normal  stress  when  the  yield  condition  (7)  holds,  one  is  led  to 

f(  (T-p  r2,  (Tz)  - V(-  tf^-pHsgnt  (T1  - r2)  + sgn(  <T^  - (T3)]/D, 

D = (^-/(-(T^pJCagnl  <Tr  r2)  + sgn  ( <7^  - (Tj)] 

(34) 

+ C2*/(”  ^”2  " P)  ^sgn  ^ ^2  " ^3)  + s6n(  I 

+ ^ '/(“  ^3  “ P)  [agn(  ^ - 0^)  + sgn  ( CT3  - 0"2)  I. 

12.  The  third  soil  theory.  In  another  possible  soil  theory 
the  yield  condition  (7)  is  used  together  with  f(  07, , (T^ ) as 
defined  in  (30) . 
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13.  The  fourth  soil  theory.  In  this  theory,  in  analogy  with 
the  third  soil  theory,  one  uses  the  yield  condition  (7)  together 

with  f(  cr1,<r2i(rz)  as  defined  ln  131). 

14.  Testing  the  theories . Reyond  this  point  experimental 
evidence  is  required  for  further  progress.  It  is  essential  to 
teat  the  various  theories  suggested  or  others.  It  is  necessary 
to  determine  the  form  of  V as  a function  of  u in  (22)  , and  the 
form  of  as  a function  of  the  plastic  work  W.  In  addition, 
the  elastic  constants  must  be  determined  for  any  given  sample. 

It  is  desirable  to  isolate  the  unknowns  in  different  testing  pro- 
grams if  possible. 

The  problem  of  determining  the  elastic  constants  can  be 
isolated.  Since 

J = ^ljdeij  " dW' 

by  (5)  , 

- flAi?  <Tn  <Tjj)  = /(7ijd€ij  - A W (35) 

The  constants  OC  and  f3  can  be  determined  if  the  stresses  and 
strains  and  the  generation  of  plastic  energy  can  be  measured  in 
the  process  of  a deformation.  The  plastic  energy  is  converted 
into  heat  and  could  presumably  be  measured  by  measuring  the  tempera- 
ture of  the  sample.  The  expressions  in  (35)  represent  lines  in 
the  , p -plane.  Experimental  errors  will  probably  prevent  a 

common  intersection  of  these  lines,  but  it  will  probably  prove 
easy  to  select  by  eye  a point  in  the  eC  , j3  - plane  suitably  near  all 
the  lines.  A more  sophisticated  least  square  or  linear  regression 
procedure  could  be  employed  if  desired. 
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The  experimental  quantities  in  equation  (35)  are  integrated 
expressions,  so  exceptional  accuracy  will  probably  not  be  required 
in  their  measurement „ However  the  determination  of  the  functional 
natures  of  V and  <f>  must  use  point-to-point  data,  so  higher 
accuracy  will  be  required  for  them.  It  may  be  desirable  to 
devise  tests  not  requiring  the  measurement  of  W,  since  this  will 
certainly  be  difficult  to  do  with  high  accuracy. 

To  avoid  complications  the  test  set-up  should  have  symmetries 
of  such  a nature  that  the  principal  directions  are  known.  Then 
by  (4),  (5),  (29), 

d(e1  - - p<r)/ f(  <rx,  <r2,  <r3)  = 

= d(e2  - occr2  - p<r)/t(  <r2,  <r3,  = (36) 

= d(e3  - <x<r3  - fl<T)/ f(  <r3,  <rlt(r2 ) = <*w 

where  (T^,  (Tg,  are  principal  stresses  and 

a~  = <rx  + <r2  + <r3.  (37) 

For  the  theories  considered  here  in  s§  9-13,  when  common  factors 
are  cancelled,  the  two  equations  on  the  left  of  (36)  have  denomina- 
tors that  are  functions  of  T*.  / T and  Ti ./  T only. 
Eliminating  the  latter  ratio  between  the  two  equations  we  may  solve 
for  . By  (7)  this  gives  a point-by  point  expression  for 

* . 

v.  being  known,  p can  be  determined  from  (1) . In 

fact 

p = (f.VT,  - <rM)/(i  - W (ss) 

gives  a point-by-point  determination  of  p. 
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15.  Static  compression  of  a cylindrical  soil  sample . We 
will  specialize  the  material  in  the  last  section  to  the  case  of 
cylindrical  soil  samples  commonly  used  in  soil  testing.  The  sample 
will  be  compressed  axially.  The  sides  may  be  unconfined,  rigidly 
confined,  or  confined  in  some  other  manner  imposing  a known 
relationship  among  the  stress  and  strain  components.  We  will 
assume  the  stress  and  strain  components  and  the  intrinsic  pressure 
to  be  constant  throughout  the  sample.  In  practice  this  will  be 
approximately  true  before  the  sample  begins  to  exude  water  except 
near  the  ends  of  the  column.  Then  by  (21),  U = const.  By  (22) 

u = q(€)  , q'(S)  < 0 (39) 

Establish  a cylindrical  f ,0  , z-coordinate  system  with 
z-axis  parallel  to  the  elements  of  the  column  ana  let  <Tr  , 0*  , G T, 
er  * ^ * ez  denote  principal  stress  and  strain  components.  Then 
Oq  = Of  = (Tp  = s <7"z  = <Tm,  e0  = er*  The  yield  conditions 
(7)  and  (14)  take  on  the  same  form.  In  the  elastic  zones,  by  (1), 

- <rr  - p > (-  <rz  - p)k,  (40) 

and  in  the  plastic  zones 

- <Tr  - P = (-  (Tz  - p)k,  (41) 

where  in  either  case 

k = (l-sin^)/(l  + sin^  ) = (1  - 77  )/(l  + 2T\)  (42) 

By  (30) - (34)  , (41), 

f(  <rz , rr,  <rr)  = -k  * /( <rr  - k1  trz) , 


(43) 
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where  5 = 1 for  the  Drucker-Prager  theory  and  the  first  soil 
theory  (which  are  identical  in  this  problem) , S = -1/2  for  the 
second  soil  theory,  and  5=0  for  the  third  and  fourth  soil 
theories  (which  are  identical  in  this  problem) . 

By  (35), 

<*Al<r?  + \<fy  + «■,))  = 

(44) 

= /(2(Trder+  J*zd62)-4w, 

from  which  0 C and  p may  be  determined. 

By  (36),  (43), 

2k*  = - d(ez  - o t<rz  + flT)/ d(er  - *<Tr  + (3T) , (45) 

from  which  k may  be  determined. 

By  (38)  or  (41), 

P = (k  <TZ  - <Tr)/(  1 - k),  (46) 

from  which  p may  be  determined. 

When  p is  constant  u may  be  determined  as  a function  of  6 
c 

up  to  an  additive  constant.  When  p varies  with  W,  samples  of  the 

c 

same  skeletal  structure  but  different  water  contents  may  be  run 
to  determine,  within  an  additive  constant,  the  function  of  W and 
the  function  of  6 which,  added  together,  equal  p, 

16,  The  case  k = const.  In  this  case  (45)  may  be  integrated, 
giving 

e((-<rr)  + pf(-<T)  + Y er  - (-e)  = const.,  (47) 

where 

p = (*y-  ^#)/(3  + y),  k*=i+|v.  (48) 
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The  constant  on  the  right  of  (47)  may  be  determined  by  evaluating 
it  at  a particular  stress-strain  state  of  the  sample.  This  may 
then  be  transposed  to  the  left  side.  Therefore  we  may  take  the 
constant  on  the  right  of  (47)  as  zero  if  instead  of  taking  the 
true  stresses  and  strains  on  the  left  of  (47),  we  take  the  differ- 
ence between  these  stresses  and  strains  and  those  of  one  particular 
stress-strain  state  of  the  sample. 

We  may  now  evaluate  Ot  , fi* , and  *V  in  (47)  by  a least 
squares  procedure.  Suppose  Yi  observations  are  made  as  the 
sample  is  compressed  and  denote  the  observed  quantities  by  adding 
a subscript  i to  the  notations  used  above,  where  i = 1,2,,..,  n. 


We  wish  to  minimize 

1 n 
v — A 


+ 13  ‘-''I1  + 


where 


= A«2  + 21  of  ^ * + 2J  eCY  + By*  2 + 2K  yS  V 
+ C y2  - 2 Fof  - 2G^S ' - 2Hy  + D, 


A = 


n 


B = i 


n 


n i=l 


(-<r,r. 


(49) 


° = H 


1 = n 


n 


J = i 


n ^ (-<rrl)erl.  (50) 


K =i 


n 


F = i 
n 


= H ^(-o-i^n- 


= H ^ (-<r1)(-e1),  h = | Z<6rl)(-V- 


n 

z 

i=l 

n 


(-  <rrl)  (-«!>, 
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of 


By  (49),  Y is 

minimized  when  oc  , fQ  , and 

V are  solutions 

A U 

+ 

ip*  + = p, 

I OC 

+ 

B^g'  + K > = G, 

(51) 

JM 

+ 

K a'  + c y = H. 

The  corresponding  value  of  Y is 

Y = D - Ftf  - Gfi'  - H y 


(52) 


By  (49)  this  quantity  may  be  judged  small  if  it  is  small  with 
respect  to  D. 

17.  An  example . The  foregoing  theory  was  applied  to  one 
case  of  a soil  sample  measurement  made  with  a Hveem  stabilometer . 
The  accuracy  of  the  test  was  below  the  requirements  of  the  theory, 
but  the  results  will  be  included  anyway  simply  as  an  illustration. 

The  values  of  the  stresses  and  strains  over  the  range  con- 
sidered are  given  in  columns  1-4  of  Table  1,  The  derivatives 
given  in  columns  7-9  were  calculated  numerically. 

After  subtracting  the  initial  values  in  columns  3-6  from  the 


other  entries,  the  quantities 
They  are 

A-H 

in  (50)  may  be  calculated 

A = 

781.000  lb2/in4, 

B = 

8750.273  lb2/in4 , 

C = 

2.951*10"6 

D = 

7.270*10"6 

I = 

2611,383  lb2/in4. 

J = 

47  .456 *10*’3  lb/in2, 

K = 

158.970e10_3lb/in2, 

P = 

75  ,067  *10”3  lb/in2, 

G = 

251 ,397  *10”3lb/in2 , 

H = 

4 .616 ’lO"6 . 
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TABLE  1. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

-^z 

"6z 

"*r 

€r 

- T 

-e 

io3de  /da; 

z'  z 

lo’deyd^  d^/d<r2 

lb/in2 

lb/in2 

lb/in2 

in2/lb 

in2/lb 

-a 

• 

03 

.0286 

21.0 

.0036 

89.8 

.0214 

.251 

.074 

.600 

55.7 

.0304 

26.0 

.0041 

107.7 

.0222 

.207 

.063 

.660 

63.7 

.0319 

31.5 

.0046 

126.7 

.0227 

.176 

.063 

.715 

71.6 

.0332 

37.5 

.0050 

146.5 

.0232 

.156 

.045 

.748 

79.6 

.0344 

43.5 

.0053 

166,6 

.0238 

.140 

.036 

.773 

87.5 

.0355 

49.5 

.0055 

186,5 

.0245 

.125 

.030 

.796 

95.5 

.0364 

56,0 

.0058 

207.5 

.0248 

.112 

.026 

.810 

103.5 

.0372 

62.5 

.0060 

228.5 

.0252 

.100 

.021 

.824 

111.4 

.0379 

69.0 

.0061 

249.5 

.0267 

.090 

.018 

.840 

The  solution  of  equations 

(51)  gives  oC 

- 18 .5  *10"6 

in2/lb , 

/»'  = 

9.35*10 

“6  in2 /lb. 

y = 

.763 . By 

(M), 

fi  = 1.26 

*10"6  in2/lb 

e 

By  (6) 

, E = 

58,  100  lb/in2, 

V = .073. 

By  (48),  k 9 = : 

1.381. 

Since  0 < k < 1,  this  indicates  the  appropriateness  of  the  second 
soil  theory  in  which  $ = - 1/2.  This  gives  k = .524,  corres- 
ponding, by  (42),  to  a shearing  resistance  angle  ^ = 18°.  By 
(52)  the  corresponding  minimum  value  of  Y is  ,036 •10~6,  so 
Y/D  = .0050. 

Using  these  values  of  cC  and  , point  to  point  values  of 
k were  computed  from  (45)  and  columns  7-9  of  Table  1.  The  results 
are  graphed  in  Pig.  2.  The  solid  line  corresponds  to  the  value 
k = .524  calculated  in  the  last  paragraph. 
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The  quantity  p may  now  be  computed  from  (46)  and  Table  1. 
The  results  are  given  in  Fig.  3 using  k = .524. 
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Appendix  . 

The  following  theorem  is  a slight  extension  of  one  due  to 
C.  B.  Morrey,  R.  M.  Lakness,  and  E.  Parzen. 

Theorem  Let  and  be  symmetric  tensors  connected 

by  a relation 

^ j = Pi j ^ Tpq^  - PiJ^l»  T22*  733*  T12>  Th*  Til* 

which  is  invariant  to  cartesian  coordinate  transformations.  Then 
there  exists  a function  f(x,y,z)  for  which 

f(x,y,z)  = f(x,z,y)  (B) 

such  that 

5lj  = *o  Sij  + Tij  + S2  Tly  T yJ  (C) 

where 

*o  = T2T3P1  + T3T1P2  + T1T2P3> 

*1  = -(Tg  + T3)Pl  - (T3  + T^Pg  - (Tx  + T2)P3,  (D) 

*2‘  P1  + P2  + P3  * 

where 

P1  = flT^Tg.TjJ/Cdj^  - TgHTj-  T3)J, 

P2  = *<WT1>/*(V  T3)(T2  - Tpj.  (E) 

P3  = f<WT2>/[<T3  - Tl)<V  Tg)  J, 

^l,^2,l^3  b«i“S  bbe  principal  values  of  the  tensor  T ^ , i.e.#  the 
roots  of  the  equation 

T3  - IXT2  + I2  T - I3  = 0 


(P) 
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where 


*1  = T*ur 


X2  ~ T11T22  + T22T33  + T33T11  " T12  ” T23 


- T 


31 


(G) 


I3  " ^T11 I . 


Proof,  A cartesian  coordinate  transformation  may  be  written 


xi  ~ aiv  yv 


where 


a*i  a«j 


~ ^lj»  ai  (t  ajfi  " ^ij 


(H) 


(I) 


The  matrices  BSijII  and  Htu«  of  tensor  components  are 

transformed  by  (H)  to  new  matrices  ||  s^  j II  and  «hjH  such 

that 


Sij  = 8rfj.Vj,l  * 


Tij  ~ ai*  aj(«  ‘ 


(J) 


Since  the  relationship  (A)  is  invariant, 


3tf( f " P*(l  ^mn^ 
Substituting  this  and  ( J)  into  (A), 


aqp  } = (tmn}  ai«aj/»  (K) 

In  particular  [8;  Gh.  I,  §3],  there  exists  a transformation, 
with  matrix  i(  ®ijll  » to  principal  axes,  i.e,,  such  that  the 
tensor  components  t^j  satisfy 

= h sij  U> 


Then  by  ( J)  , (K)  , taking  = a^, 

Tij  = T>  ai v aj>  * 

Pij^Tpq^  = f *<•  (T1»T2»T3)  ai*  aj^  * 


where 


j^Tl,T2,T3^  = Pij^Tm  & mn^ 

Now  let  t^j  = t^j  in  (K) , 

Fi ap*aq*^  = f *(»  ^Tl,T2»T3^ai«C  aj|f 


Insert 


iiau/l  - 


Then 


0 10 
0 0 1 
1 0 01 


|T2  o 0 


HT*ai*ajJl  = ||°  T3  0 

0 0 T, 


Therefore  by  (0)  and  (P), 


f22 

f23 

f12 

Hflj(T8'T3-Tl)H  = 

f23 

f33 

f13 

f12 

f13 

fll 

(M) 

(N) 

(O) 
(P) 


the  arguments  of  the  functions  f^  on  the  right  being  t1,tg,t3 
This  relation  gives 


Now  insert 


in  (P).  This  gives 


t a.  a 
a ia 


so  by  (0)  and  (P) , 


||  fiJ(tl»t3,t2)| 


10  0 
0 0 1 
0-10 

t 0 0 

0 t^  0 


1 

0 

0 

*z| 

fll 

f13 

"fl2 

mm 

f13 

f33 

”f23 

"f12 

“f23 

f22 

This  gives 

fll^tl,t2,t3^  = fll*tl,fc3,t2^  f12^t1,t2,t3^  = "fi3^ti,t3,t2^ 

(R) 

f13(tl»t2»  fe3)  * **12  ^1*^3*  ^2^  • 


Defining 


**U,y,z)  = f11(x,y,z). 


(B)  follows  from  (R).  Also,  from  (R),  f12  = 0.  Therefore  by  (Q), 

fi3  = °*  f23  ■ °*  and  f22^x,y,z^  “ f33^x»y»z)  s f(z»*»y)« 

To  complete  the  proof  we  must  show  that  the  right-hand  sides 
of  (C)  and  (N)  are  equa  1.  Inserting  (M)  into  (C),  and  using  (I), 
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Now  by  (D)  , 

#0  + *1T1  + *2T12  = <V  T2>‘T1  - T3>P1  : 

*0  + #1T2  + *2T22  = <T2  * T3>(T2  * T1>P2 
#0  + #^3  + *2T32  = (I,  - TjHTj  - Tg)P3 


« W3> 

= f(Tg,T3,T1) 

= ^WT2> 


Substituting  this  into  (S) , we  obtain  the  right-hand  side  of 
(N) , completing  the  proof. 
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